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Abstract 

We applied an improved Coulomb correction method developed by us recently to data on identical 
AA-pairs production in S + Pb and p + Pb reactions at 200 GeV/c obtained by NA44 Collaboration. 
To analyse the whole range of the momentum transfers measured the method of ’’seamless htting” 
has been proposed and used together with the asymptotic expansion formula for the Coulomb wave 
function. We found that such Coulomb corrections lead sometimes to different than previously re¬ 
ported (by NA44 Collaboration) interaction region and strongly influence the long range correlations. 
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Introduction: Recently NA44 Collaboration has reported their data on the Bose-Einstein correlations 
(BEC) of pairs produced in S -I- Pb and p -I- Pb reactions at 200 GeV/c j|]. In our previous 

work we have analysed these data by making use of the various source functions with the long range 
correlation (without, however, invoking any sort of Coulomb corrections). The data Q have been cor¬ 
rected for Coulomb interactions by applying only the Gamow factor. As was pointed out some time ago 
by Bowler this is, however, not sufficient |^. In our recent works il the still improved method of 
Coulomb corrections was presented but not yet applied to any concrete data set. In the present Letter, 
we apply it therefore to the analysis of data for pairs production in S -I- Pb and p -j- Pb reactions 
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mentioned above (in the whole measured momentum transfer region) and compare our results with those 
obtained before in |^ . To be able to analyse the whole range of measured momentum transfer Q and to 
avoid wild oscillations developing at large Q’s (cf. Fig. la) we have to use the the asymptotic expansion 
of the Coulomb wave function together with the procedure of ’’seamless fitting” (SF) explained below 
(cf. Figs, lb and 2). 


Theoretical formula of BEC with Coulomb wave function: To write down an amplitude A 12 
satisfying Bose-Einstein statistics it is convenient to decompose the wave function of identical (charged in 
our case) bosons with momenta pi and p 2 into the wave function of the center-of-mass system (c.m.) with 
total momentum P = ^{pi +P 2 ) and the inner wave function with relative momentum Q = {pi —P 2 ) = 2(?. 
It allows us to express A 12 in terms of the confluent hypergeometric function $ [Q: 

^12 = ;^[^'(q,r)-f 4 's(q,r)] , ( 1 ) 

dt(q, r) = r(l -I- i?7)e“'"'''/^e*‘* '’<i>(—1?7; 1; i(7r(l — COS 0 )) , 

'I's(q, r) = + + 


where r = xi — X 2 and the parameter rj = ma/2q. Assuming factorization in the source functions, 
p{i"i,r 2 ) = p(ri)p(r 2 ) = p{R)p{r) (here R = ^{xi + X 2 )), one obtains the following expression for 
theoretical BEC formula including the improved Coulomb correction ||]: 
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where G{q) = 2Trp/(e^'^'^ — 1) denotes Gamow factor and the first and the second parentheses in eq. (i 
correspond to the first and the second terms in eq. (||)0 and 

r(zr7-|-n) 1 
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To analyse data corrected only by the Gamow factor using our formulae we should use the following 
ratio: 


^(±±:GC)/^BG(q ^ 2q) = Rcc/G{q) 


= c(l -I- Alc -I- Aec) 


1 + A 


F^2B 


1 -I- Alc -I- Aec 


(l + pQ). 


^For the exact formulae for Ag(;;; and F^2B 0* 


(4) 
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It should be noted that the normalization and an effective degree of coherence, i.e., the denominator of 
the ratio -B 2 b /(1 + Aic + Aec), are related to each other. Notice also that other parameters like the 
additional normalization factor c, the long range correlation 7 and A are introduced by hand. 


Source function: To obtain an explicit expression, we have to decide on some form of the source 
function. In the present Letter, we shall use the Gaussian source distribution, p{r) = exp(— 

For this type of source function we have the following formulae for the elements of eqs. and i) i: 
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E 2 B = exp ( - — 


To analyze data corrected by the Coulomb wave function as was done in 
formula (Q) replacing it by the following one: 


(5) 

( 6 ) 

we should modify the 


^(±±:CC)/^BG(g ^ 2g) = 


Rcc 


^(9)(l + Aic + Aec) 

-E'2B 


1 + A- 


(1 + 7<3)- 


(7) 


I + Aic + Aec 

For the sake of reference, we write down here also the conventional formula (i.e., a kind of standard 
formula without corrections due to the final state interactions): 

^(±±:Standard) ^ 2q) = c[l + XE2 b] (1 + iQ)- ( 8 ) 


Asymptotic expansion of the Coulomb wave function: First of all, it should be noted that 
the expansion in eq. (^) has to be the limited to giimit(= Q/ 2 ) only due to mathematical properties of the 
confluent hypergeometric function used [Q|. This can be seen as wild oscillation developing in Fig. 1 (a) 
where the eq. (||) has been simply used. If we, instead, set the Coulomb correction to zero in the region 
q > ®imit limit, a small step appears as seen in Fig. I (b). Therefore in order to analyse in a consistent way 
the whole region of the momentum transfer measured we have to use the following asymptotic expansion 
of the Coulomb wave function: 


^ [d/-^“(q, r) + r)] 

'k’^’’ym(q, r) = exp{i((7z — 77ln(r — z))} X (i + t 


iq{r - z) 


+f{0) 


exp{*(qr — r]ln{2qr))} 


'I'g’’y™(q, r) = exp{i{—qz — rj ln(r + z))} x ( 1 + t 


iq{r + z) 


+/(7r- 6») 


exp{*(qr — p ln( 2 < 7 r))} 


(9) 


( 10 ) 


( 11 ) 


where z = r cos 9 and 


/(^) = ®^P{“ 2 ir 7 lnsin( 6 l/ 2 ) + 2i argr(l + irj)}. 

sm ( 6 // 2 ) 
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In analyses we should assure a smooth connection between both regions of g = Q/2. To avoid the diver¬ 
gence of denominators (1 ± cosd), we have to introduce a cutoff parameter e (of the order of e ~ 10 “^) 
such that (1 ± COS0) > e always^. This procedure, shown as a flow chart in fig. 2, is called the “seamless 
fitting (SF)”. 

Analyses of data by SF method: Our results obtained in terms of the new formula (Q) are shown 
in Table I and Figs. 3 and 4. Whereas the parameter 7 (representing here the influence of long range 
correlations) increases now noticeably in comparison with that obtained previously in we find that 
the interaction region represented by i? = 1/2/3 remains S -I- Pb reactions almost the same. Only in p 
+ Pb collisions the estimated R becomes larger with inclusion of Coulomb correction than that obtained 
in |Q|. These facts suggest that we should be careful in interpreting any data (at least those for kaons) 
which were corrected for Coulomb interactions only by the Gamow factor. 

Concluding remarks: We have proposed the possible method of applying the Coulomb correction 
for the BEC in the whole region of momentum transfer, the “seamless fitting” (SF). We confirm that this 
method works well when analysing data corrected only by the Gamow factor. 

Our analyses of data of KK pairs in S -I- Pb reaction Ql shows (cf. Table I) that R{K'^K~^) = 4 fm 
and R{K~K~) = 3 fm (i.e., they differ substantially), contrary to the estimation provided in [|| that 
R{K^K^) « 3 fm. This is an important result for the study of signals of the Quark Gluon Plasma 
(QGP) (see refs. |^, ^, ^). Moreover, we found that the long range correlations are strongly affected by 
Coulomb corrections (most probably because of the long range character of Coulomb interactions). 

For the sake of reference, we show also in Table II and Fig. 5 results of our analysis of data for S -I- 
Pb ^ TTTT -I- X reaction Q (which were corrected by the Coulomb wave function method |^) performed 
by using both eq. (^ and the “standard” formula (eq. (^). As one can see there are no significant 
differences between parameters estimated by means of these two formulae, in particular the magnitude 
of the interaction region is in both cases almost the same. 

Acknowledgements: The authors would like to thank S. Esumi, T. Nishimura, and S. D. Pandey 
(members of NA44 Collaboration) for useful conversations and correspondences. This work is partially 
supported by Japanese Grant-in-Aid for Scientihc Research from the Ministry of Education, Science and 
Culture (# 06640383). 

^We have confirmed that the parameter e depends on the magnitude of the interaction region (R). 

®This shows that Coulomb corrections can be important in the bigger range of momentum transfers than considered so 
far and justifies a posteriori our present investigation. 
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Table I: Parameters of kaonic BEG obtained from S -I- Pb and p + Pb reactions. 


reaction 

formula 

c 

R [fm] 

A 

7 

xVndf 


standard 

0.998±0.008 

3.384±0.171 

1.010±0.071 

— 

58.0/33 

S -b Pb ^ 


1.101±0.022 

3.991±0.234 

0.977±0.078 

-0.420±0.079 

36.8/32 

K+K+ + X 

eq. (4) 

0.966±0.007 

0.987±0.015 

3.918±0.000 

4.285±0.000 

0.804±0.059 

0.825±0.074 

-0.071±0.077 

30.9/33 

30.9/32 


standard 

0.988±0.021 

3.058±0.285 

1.070±0.142 

— 

31.9/32 

S + Pb ^ 


0.989±0.066 

3.061±0.345 

1.069±0.159 

-0.004±0.279 

31.9/31 

K-R- PX 

eq. (4) 

0.944±0.017 

0.872±0.048 

3.089±0.000 

3.066±0.012 

0.836±0.111 

1.030T0.178 

0.411±0.277 

35.0/32 

32.5/31 


standard 

0.989±0.027 

2.134±0.291 

0.636±0.089 

— 

34.3/26 

p -b Pb ^ 


1.149±0.078 

2.698±0.559 

0.493±0.114 

-0.546±0.237 

30.9/25 

K+K+ + X 

eq. (4) 

0.975±0.019 

0.996±0.061 

2.731T0.000 

2.310±0.001 

0.520±0.091 

0.431±0.133 

-0.180±0.256 

32.7/26 

30.7/25 


Figure Captions 

Fig. 1. (a) Analyses of data for S -f Pb ^ K~K~ + X reaction by eq. (^. (b) The same but with 

Coulomb corrections switched off for q > = Ql‘^ limit. (Results for pair production 
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Table II: Parameters of S + Pb ^ 7r=*=7r=*= + X reactions. 


formula 

c 

R [fm] 

A 

7 

xVndf 

standard 

0.800±0.003 

0.824±0.010 

4.506±0.327 

5.0I5±0.421 

0.461±0.042 

0.450±0.044 

-0.178±0.065 

17.5/16 

10.9/15 

eq. (7) 

0.800±0.003 

0.825±0.010 

4.575±0.335 

5.105±0.436 

0.509±0.048 

0.502±0.052 

-0.180±0.065 

17.8/16 

11.0/15 


looks similar with only change being in the value of gumit which depend on the value of radius 
parameter R.) 

Fig. 2. Flow chart for our procedure of “seamless fitting (SF)”. 

Fig. 3. Results of SF for BEC for kaons produced in S + Pb collisions; (a) for pairs; (b) for 

K~K~ pairs. 

Fig. 4. Results of SF for BEC for kaons produced in p + Pb ^ K K + X reaction. 

Fig. 5. Results of SF for BEC for pions produced in S + Pb 7r+7r+ + X reaction. 
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This figure "fig 1-1.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9507293v2 



This figure "fig2-l.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9507293v2 



This figure "fig 1-2.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9507293v2 



